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How to use this workbook 


This workbook is designed to be used in conjunction with the authors free online video tutorials. 
Inside this workbook each chapter is divided into learning modules (subsections), each having its own 
dedicated video tutorial. 

View the online video via the hyperlink located at the top of the page of each learning module, with 
workbook and paper / tablet at the ready. Or click on the Engineering Mathematics YouTube Workbook 
playlist where all the videos for the workbook are located in chronological order: 

Engineering Mathematics YouTube Workbook Playlist 
http://www.youtube.com/playlist?list=PL13760D87FA88691D [YTPL]. 

While watching each video, fill in the spaces provided after each example in the workbook and annotate 
to the associated text. You can also access the above via the authors YouTube channel 

Dr Chris Tisdells YouTube Channel 
http://www.youtube.com/DrChrisTisdell 

The delivery method for each learning module in the workbook is as follows: 

• Briefly motivate the topic under consideration; 

• Carefully discuss a concrete example; 

• Mention how the ideas generalize; 

• Provide a few exercises (with answers) for the reader to try. 

Incorporating YouTube as an educational tool means enhanced eLearning benefits, for example, the 
student can easily control the delivery of learning by pausing, rewinding (or fast-forwarding) the video 
as needed. 

The subject material is based on the authors lectures to engineering students at UNSW, Sydney. The 
style is informal. It is anticipated that most readers will use this workbook as a revision tool and have 
their own set of problems to solve - this is one reason why the number of exercises herein are limited. 

Two semesters of calculus is an essential prerequisite for anyone using this workbook. 
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About the author 

“With more than a million YouTube hits, Dr Chris Tisdell is the equivalent of a best-selling author or 
chart-topping musician. And the unlikely subject of this mass popularity? University mathematics” 
[Sydney Morning Herald, 14/6/2012 http://ow.ly/o7gti l . 

Chris Tisdell has been inspiring, motivating and engaging large mathematics classes at UNSW, Sydney 
for over a decade. His lectures are performance-like, with emphasis on contextualisation, clarity in 
presentation and a strong connection between student and teacher. 

He blends the live experience with out-of-class learning, underpinned by flexibility, sharing and openness. 
Enabling this has been his creation, freely sharing and management of future-oriented online learning 
resources, known as Open Educational Resources (OER). They are designed to empower learners by 
granting them unlimited access to knowledge at a time, location and pace that suits their needs. This 
includes: hundreds of YouTube educational videos of his lectures and tutorials; an etextbook with each 
section strategically linked with his online videos; and live interactive classes streamed over the internet. 

His approach has changed the way students learn mathematics, moving from a traditional closed 
classroom environment to an open, flexible and forward-looking learning model. 

Indicators of esteem include: a prestigious educational partnership with Google; an etextbook with 
over 500,000 unique downloads; mathematics videos enjoying millions of hits from over 200 countries; 
a UNSW Vice-Chancellors Award for Teaching Excellence; and 100% student satisfaction rating in 
teaching surveys across 15 different courses at UNSW over eight years. 

Chris has been an educational consultant to The Australian Broadcasting Corporation and has advised 
the Chief Scientist of Australia on educational policy. 

Subscribe to his YouTube channel for more http://www.youtube.com/DrChrisTisdell . 
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1 Partial derivatives & applications 

1.1 Partial derivatives & partial differential equations 

View this lesson at http://www.youtube.com/watch?v=SV56UC31rbk , [PDE]. 



Motivation 





Partial differential equations (PDEs) are very important in modelling as their solutions unlock the 
secrets to a range of important phenomena in engineering and physics. The PDE known as the wave 
equation models sound waves, light waves and water waves. It arises in fields such as acoustics, 
^electromagnetics and fluid dynamics. 


Example 


Consider the wave equation 

d 2 w 2 d 2 w 
= cr 


c is a constant. 


dt 2 dx 2 5 

Show that a solution is given by 

w(t , x) := cos(2x + 2 ct). 



The bigger picture 





• When verifying solutions to PDEs, compute the necessary derivatives and show that the 
PDE holds. 

• In certain circumstances, the chain rule will be required to compute the partial 
derivatives of interest. 
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Exercises 


< 



Consider the wave equation 

d 2 w 2 d 2 w 

tt- = c ——, c is a constant. 

ot z ox 1 

Show that another solution is given by 

w(t , x ) = ln(3x + 3 ct). 

Show that the wave equation has solutions of a more general form 

w(t , x) = f{x + ct) 

^where/is a sufficiently smooth function. (Hint: use the chain rule.) 


1.2 Partial derivatives & chain rule 

View this lesson at http://www.youtube.com/watch?v=HOzMR22HsiA , [Chain]. 



Motivation 





The chain rule is an important technique for computing derivatives and better-understanding rates 
^of change of functions. For functions of two (or more) variables, the chain rule takes various forms. ^ 



Example 


s' 



Let the function w ._ ^ have continuous partial derivatives. If we make a change of variables: 

x — r cos Q ; y = r sin 6 then show that 


V. 


dw 

dr 


d f a ,df ■ n 

—— cos u + —— sin u. 
ox oy 


J 


r 


The bigger picture 


"V 


• It is generally helpful to draw a simple diagram to understand which form of the chain 
rule to apply. 

• For functions of two (or more) variables, the chain rule takes a more profound form than 
the one-variable case. 

v___ J 
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Exercises 


< 



Let the function w f(x,y) have continuous partial derivatives. If we make a change of variables: 
x — r cos 6 ; y — r sin 6 then show that 


V 


dw 

~d6 


df . a df 

— —r sm 0 + —— r cos u. 
ox oy 




1.3 Taylor polynomial approximations: two variables 

View this lesson at http://www.youtube.com/watch?v=ez HZZ9H2ao , [Tay]. 


s — Motivation - n 

Taylor polynomials are a very simple and useful way of approximating complicated functions. Taylor 
polynomials are desirable types of approximations as their polynomial structure make them easy to 
work with. 

V_ J 
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Example 




A 

Using an appropriate Taylor polynomial, compute an approximation to 


\ 

\J (1.02) 3 + (1.97) 3 . 

J 


The bigger picture 


The Taylor polynomial up to and including linear terms is usually a good starting 
approximation (unless a greater degree of accuracy is required) 

T(x,y) := f(a,b ) + f x (a, b)(x - a) + f y (a,b)(y - b). 

The use of a linear Taylor polynomial approximation geometrically equates to 
approximating surfaces by a suitable tangent plane. 



Exercises 





Using an appropriate Taylor polynomial, compute an approximation to 


a/ (0.98) 3 + (2.03) 3 . 


[Ans: 2.9]. 

V___ 

1.4 Error estimation 

View this lesson at http://www.youtube.com/watch?v=wLrHCkesOA8 , [Est]. 



Motivation 





When taking measurements (say, some physical dimensions), errors in the recorded measurements are 
a fact of life. In many cases we require measurements to be within some prescribed degree of accuracy. 

We now look at the effects of small variations in measurements and estimate the errors involved. 

V___ 



Example 





A cylindrical can has height h and radius r. We measure the height and radius and obtain 12cm and 
5cm respectively, with errors in our measurements being no more than 0.5mm. As a result of the 
errors in our measurements for h and r, obtain an estimate on the percentage error in calculating 
the volume of the can. 

V_ J 
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r 


The bigger picture 


• The main inequality used can be derived from the linear Taylor polynomial 
approximation for the function involved. 

• Don’t forget to unify the units involved! 




J 



Exercises 





• Re-examine the previous example by switching the recorded measurements for height 

and radius and then estimate the percentage error in V. Compare your estimate with the 
estimate of the previous example. [Ans 11/6%] 

• We measure the length l and width w of a plot of land to be 40m and 10m respectively, 
with the errors in these measurements being no more than 0.5cm. As a result of the 
errors in l and w, obtain an estimate on the percentage error in calculating the area of the 
plot.[Ans: 1/16%] ^ 

1.5 Differentiate under integral signs: Leibniz rule 

View this lesson at http://www.youtube.com/watch?v=cipIYpDu9YY , [Leib]. 



Motivation 



\ 


Leibniz rule: 

I(x) := f f(t,x) dt 

J a 

has derivative 

I'(x) — [ -J-(t,x) dt. 

Ja 9x 

There are two important motivational points concerning Leibniz’ rule: 

• If a function I defined by an integral is useful for modelling purposes, then it would seem 
to make sense that the derivative T would also give us insight into the problem under 
consideration; 

• Leibniz’ rule can be applied to evaluate very challenging integrals. 
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Some max/min problems for multivariable functions 


2 Some max/min problems for 
multivariable functions 


2.1 How to determine & classify critical points 

View this lesson at http://www.youtube.com/watch?v=LGFlLcLnSfE , [CPI]. 



Motivation 





Functions of two (or more) variables enable us to model complicated phenomena in a more accurate 
way than, for example, by using functions of one variable. The determination and classification of 
critical points of functions is very important in engineering and the applied sciences as this information 
^is sought in a wide range of problems involving modelling. 



Example 





Determine and classify all the critical points of the function 

f(x, y) := x 3 + y 3 — 3 xy + 15. 

v_ J 
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The bigger picture 


Assume that/is continuous and has continuous partial derivatives up to and including all those of 
the second-order. 

• Determine the critical points of/by solving the simultaneous equations :f x = 0 and^ = 0. 

• If (a, b) is a critical point determined from above and defining 

D(a, b ) := f xx (a, b)f m (a, b ) - [f xy (a, b)f 

then the nature of the critical points can be then determined via the 2nd-derivative test: 

1) If D(a, b) > 0 and f xx (a, b) < 0 then/has a local maximum at (a, b ); 

2) If D(a, b) > 0 and f xx (a, b) > 0 then/has a local minimum at ( a , b ); 

3) If D(a, b) < 0 then/has a saddle point at (a, b ); 

4) If D(a y b) = 0 then the 2nd-derivative test cannot be used and some other method or 

technique must be sought to classify our critical points. 



Exercises 

v_ ) 





Determine and classify all the critical points of the function 


/(/, t/) :— x ^ xy + 3x *2y -\- 5. 

^[Ans: This/has one critical point at (-2, 1) with our/having a saddle point at (-2, 1).] 


2.2 More on determining & classifying critical points 

View this lesson at http://www.youtube.com/watch?v=lq2kdnRY5h8 , [CP2]. 



Motivation 





The determination and classification of critical points of functions is very important in engineering 
and the applied sciences as this information is sought in a wide range of problems involving modelling. 
^We now look at a more involved example. 
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Example 


< 



Let A > 0 be a constant. Determine and classify all the critical points of the function 


/0, y) ■■= X 2 - Axy + y 2 + 7. 

^Hint - Discuss the cases: A = 2; 0 < A < 2; A > 2. 



The bigger picture 





Assume that/is continuous and has continuous partial derivatives up to and including all those of 
the second-order. 

• Determine the critical points of/by solving the simultaneous equations :f x = 0 and^ = 0. 

• If ( a , b) is a critical point determined from above and defining 

D(a, b ) := f xx (a, b)f yy (a, b) - [f xy (a, b)f 

then the nature of the critical points can be then determined via the 2nd-derivative test: 

1) If D(a, b) > 0 and f xx (a, b) < 0 then/has a local maximum at (a, b); 

2) If D(a, b) > 0 and f xx (a, b) > 0 then/has a local minimum at ( a , b); 

3) If D(a y b) < 0 then/has a saddle point at (a, b); 

4) If D(a, b) = 0 then the 2nd-derivative test cannot be used and some other method or 
technique must be sought to classify our critical points. 



Exercises 




Determine and classify all the critical points of the functions: 


(a) f(x, 

(b) f{x, 

y) := x 2 + xy + y 2 + 3y + 3; 

V) :=-x 2 -l-y 2 . 


[Ans: (a) This/has one critical point at (1, -2) with/having a local minimum at (1, -2). 
(b) This/has one critical point at (0, 0) with/having a local maximum at (0, 0).] 
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2.3 The method of Lagrange multipliers 

View this lesson at http://www.youtube.com/watch?v=4E-uLaRhrcA , [Lagl]. 



Motivation 


< 



The method of Lagrange multipliers is a very powerful technique enabling us to maximize or minimize 
a function that is subject to a constraint. Such kinds of problems frequently arise in engineering and 
applied mathematics, eg, designing a cylindrical silo to maximize its volume subject to a certain fixed 
^amount of building material. 



Example 





Consider a thin, metal plate that occupies the region in the XT-plane 


tt := {(x, y) : x 2 + y 2 < 25}. 

If f(x, y ) := 4.x 2 — 4 xy + y 2 denotes the temperature (in degrees C) at any point (x, y) in D, then 
^determine the highest and lowest temperatures on the edge of the plate. 
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The bigger picture 


Suppose/andghave continuous partial derivatives and: dg/dx ^ 0; dg/dy =/ 0; wheng= 0. 

• If/is our function to maximize/minimize and g = 0 is our constraint then the method of 
Lagrange multipliers involves calculating the critical points of the Lagrangian function 

L := f - Xg 

where X is the Lagrange multiplier. 

• To solve the resulting simultaneous equations, some creativity is often required. 

• If L has more than one crticial point then compare values of/at each critical point to 
determine which gives max/min values off 



Exercises 


/ 



Determine the maximum value of f(x, y) := xy subject to the constraint g(x, y) := x + y - 16 = 0. 

[Ans: At (8, 8),/has a maximum value of 64.] ^ 


2.4 Another example on Lagrange multipliers 

View this lesson at http://www.youtube.com/watch?v=5w-blyU9hy4 , [Lag2]. 



Motivation 





The method of Lagrange multipliers is a very powerful technique enabling us to maximize or minimize 
a function that is subject to a constraint. Such kinds of problems frequently arise in engineering and 
^applied mathematics. 



Example 

v_ J 





The temperature at a point (x, y) on a metal plate in the XT-plane is given by 


T(x, y) := 6xy. 

JJse the method of Lagrange to determine the maximum temperature on the circle x 2 + y 2 = 8. 
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The bigger picture 


Suppose/andghave continuous partial derivatives and: dg/dx ^ 0; dg/dy =/ 0; wheng= 0. 

• If/is our function to maximize/minimize and g = 0 is our constraint then the method of 
Lagrange multipliers involves calculating the critical points of the Lagrangian function 

L := f - Xg 

where X is the Lagrange multiplier. 

• To solve the resulting simultaneous equations, some creativity is often required. 

• If L has more than one crticial point then compare values of/at each critical point to 
determine which gives max/min values off 


Exercises 


Minimize /(x, y ) z) := x 2 + y 2 + z 2 subject to the constraint g{pc, y) := 2x — y + z — 1 = 0. 

[Ans: At (1/3, 1/6, -1/6),/has a maximum value of 2/9.] 


2.5 More on Lagrange multipliers: 2 constraints 

View this lesson at http://www.youtube.com/watch?v=7fRYd8PKeGY , [Lag3]. 



Motivation 





The method of Lagrange multipliers is a very powerful technique enabling us to maximize or minimize 
a function that is subject to a constraint. Sometimes we have two (or more) constraints. Such kinds 
^of problems frequently arise in engineering and applied mathematics. 



Example 





Maximize the function /(x, y , z) := x 2 + 2 y — z 2 subject to the constraints: 

9i(x, y , z) := 2x - y = 0; g 2 (x , y,z):=y + z = 0. 

V___ J 
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The bigger picture 


• If/is our function to maximize/minimize and g = 0, g 2 = 0 are our constraints then the 
method of Lagrange multipliers involves calculating and comparing the critical points of 
the Lagrangian function 

L f ~ [Aigi + A 2 // 2 ] 

where X , X 2 are the Lagrange multipliers. 

• To solve the resulting simultaneous equations, some creativity is often required. 

• If L has more than one crticial point then compare values of/at each critical point to 
determine which gives max/min values off 



Exercises 


< 



Determine the minimum value of / (x, y, z) := x 2 + y 2 + z 2 subject to the constraints: 

gi(x, y,z) := x + 2y + 3z — 6 = 0; ^ 2(^5 z) := x + 3y + 9z — 9 = 0. 


[Ans: Min value occurs at (1 + 22/59, 2 + 5/59, 9/59).] 
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3 A glimpse at vector calculus 

3.1 Vector functions of one variable 

View this lesson at http://www.youtube.com/watch?v=hmNlDukQ8Yk , [VecFun]. 


Motivation 


Vector-valued functions allow us more flexibility in the modelling of phenomena in two and three 
dimensions, such as the orbits of planets. The most basic of vector-valued functions are those involving 
one variable which can be used to describe and analyze curves in space. 
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Example 


< 



Consider the function r(£) := (cost, sin£, t) for t >0. 

1. Sketch and describe the curve associated with r. 

2. If a particle travels along this curve (with t representing time) then calculate r’(t) and 
show that its speed is constant. 

3. Calculate 

r(t) dt. 

\ _ ) 



r 


The bigger picture 




• Differentiation, integration and limits of vector-valued functions of one variable involve 
differentiation, integration and limits of each component function. 

• For smooth curves, r \t) is a tangent vector to the curve and points in the direction of 
motion r'(t) / \\r'(t) || . 



Exercises 





Consider the function r(t) := (2 cos t, 3 sin t) for t >0. 

1. Sketch and describe the curve associated with r. 

2. If a particle travels along this curve (with t representing time) then calculate r \t). Is the 
speed constant? 

3. Calculate 

/ 7T 

r (t) dt. 

-7T 

[Ans: 1. The curve is the ellipse x 2 /4 + y 2 /9 = 1; 2. r'(£) = (—2sin£, 3 cos t). No.; 3. (0, 0).] 
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3.2 The gradient field of a function 

View this lesson at http://www.youtube.com/watch?v=Xw5wsWDtlCg , [Grad]. 



Motivation 





The gradient field of a function is one of the basic concepts of vector calculus. It can be used to 
construct normal vectors to curves and surfaces, it can be applied to calculate slopes of tangent lines 
to surfaces in any direction, and can be very useful when integrating over curves (line integrals) to 
^compute work done. 


Example 


Consider the function f(x,y) := x * 1 2 3 + 4 y 2 . 

1. Compute V/. 

2. Show that V/is normal to the level curve /(x, y) = 16. 

3. Calculate the rate of change (directional derivative) of/at (1, 1) in the direction 
u — (1, 1). 



The bigger picture 




• The gradient readily generalizes to equations involving three or more variables, 

V/ = (. fxjyjz )• 

• For differentiable functions/the vector V/ will always be normal to level curves / level 
surfaces off 

• For differentiable functions of two variables the directional derivative of/at P(x o, yf) in 
the direction of u is just the slope of the tangent line to the surface of/with the tangent 
line lying in the vertical plane that contains (xq, Ho, /(xq, yo)) and u. 

• Directional derivatives are generalizations of partial derivatives. 



Exercises 





Consider the function f{x,y,z) x 2 + y 2 — z 2 . 

1. Compute V/ 

2. Compute a normal vector to the surface/(x, y, z) = -7 at (1, 1, 3). 

3. Calculate the rate of change (directional derivative) of/at (1, 1, 1) in the direction 
u = (1, 1, 1). 

[Ans: 1. (2x, 2y -2 z); 2 (2, 2, -6); 3. 2/V3.] 
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3.3 The divergence of a vector field 

View this lesson at http://www.voutube.com/watch?v=quZlfp59iIg , [Div]. 



Motivation 


< 



The divergence of a vector field is one of the basic concepts of vector calculus. It measures expansion 
and compression of a vector field and can be very useful when integrating over curves (line integrals) 
^and surfaces to compute flux. 



Example 


< 



a) If F(x, y, z) := ( x 2 — y, y + z, z 2 — x) then compute V*F at (1, 2, 3). 

b) If G(x, y) (0, x) then compute V • G. Sketch G and show that there is a zero net 
outflow over each rectangle. 



A cate-Lucent 


www.alcatel-lucent.com/careers 


What if 
you could 
build your 
future and 
create the 
future? 




One generation’s transformation is the next’s status quo. 
In the near future, people may soon think it’s strange that 
devices ever had to be “plugged in.” To obtain that status, there 

needs to be “The Shift”. 


27 

Download free eBooks at bookboon.com 




























Engineering Mathematics: YouTube Workbook 


A glimpse at vector calculus 


The bigger picture 


• The divergence measures net outflow of a vector field. 

• If the divergence is positive everywhere, then there is a net outflow over every closed 
curve / surface. 

• If the divergence is negative everywhere, then there is a net inflow over every closed 
curve / surface. 

• A vector field with zero divergence everywhere is called “incompressible” with zero net 
outflow over every closed curve / surface. 



Exercises 


< 



Consider F(x, y, z) := (xcosy, y sinx, z 2 — 1) and G(x, y) := (x, y ). 

1. Compute V • F and V • G. 

2. Sketch G and show that there is a net outflow over each circle in the plane. 

[Ans: 1. cos y + sin x + 2z , 2; 2. Sketch G and graphically show the net outflow is positive.] ^ 

3.4 The curl of a vector field 

View this lesson at http://www.youtube.com/watch?v=qeOTOYufpSY , [Curl]. 



Motivation 



A 


The curl of a vector field is one of the basic concepts of vector calculus. It measures rotation in a vector 
field and can be very useful when integrating over curves (line integrals) and surfaces. 



Example 





a) If F(x, y, z) := (x 2 — y,y + z, z 2 — x) then compute V x F at (1, 2, 3). 

b) If G(x, y) := (y, 0) then compute the scalar curl of G. Sketch G and discuss the 
rotation in G. 
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r 


The bigger picture 




• The curl measures rotation in a vector field. 

• If the scalar curl is positive everywhere in the plane, then there is an anticlockwise 
rotation of a paddlewheel in the plane about its axis. 

• If the scalar curl is negative everywhere in the plane, then there is clockwise rotation of a 
paddlewheel in the plane about its axis. 

• A vector field with zero curl everywhere is called “irrotational”. In the plane this means a 
paddlewheel in the vector field not spinning about its axis. 



Exercises 


< 



Consider F(x, y, z) := (zcosy, z sinx, z 2 — y 3 ) and G(x, y) := (x, y). 

1. Compute V x F and the scalar curl of G. 

2. Sketch G and graphically show that G is irrotational. 


[Ans: 1. (—3 y 2 — sin x, cos y , z cos x + z sin y ), 0; 

^_ 2. Sketch G and graphically show a paddlewheel cannot rotate about its axis.] ^ 

3.5 Introduction to line integrals 

View this lesson at http://www.youtube.com/watch?v=6YuyBnXBFXg , [linel]. 



Motivation 


s' 



A line integral involves the integration of functions over curves. Applications include: calculating work 
^done; determining the total mass and center of mass of thin wires; and also finding flux over curves. ^ 



Example 

V_ J 


< 



a) If /(x, y) := x + y and C is the unit circle with center at (0, 0) then compute the line 
integral 

J f ds. 

b) If F(x, y) := {—y + l)i + xj and C is oriented via an anticlockwise rotation then 
compute the line integral 

F • dr. 

\ _ ) 
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— The bigger picture 


For a smooth curve C with parametrization r = r (t), a < t < b 

f fds= [ b f(rmwmdt . 

JC Ja 

For an oriented smooth curve C with parametrization r = r (t), a < t < b 


/ F • dr — / F • T ds 

Jc Jc 


where T is the unit tangent vector in the direction of motion. 





In the past four years we have drilled 

* 


81,000 km 

A 


That's more than twice around the world. 



Whn am wp? fHSHHHH 


P 

We are the world's leading oilfield services company. Working 1 


globally—often in remote and challenging locations—we invent, 
design, engineer, manufacture, apply, and maintain technology 
to help customers find and produce oil and gas safely. 



Who are we looking for? 

We offer countless opportunities in the following domains: 

■ Engineering, Research, and Operations ^ 

■ Geoscience and Petrotechnical 

■ Commercial and Business 

A ^ 


If you are a self-motivated graduate looking for a dynamic career, 
apply to join our team. 

What will you be? 

careers.slb.com 

Schlumberger 
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Exercises 


< 



Let C be part of a parabola y = x 2 parametrized by r(t) = ti + t 2 j, 0 < t < 2. 

3. If /(#, y) := x 2 — y + x then compute 

J f ds. 

4. If C is oriented in the direction from (0, 0) to (1, 1) and if F(x, y) := 2xi + yj then 
compute 


L 


F • dr. 

What answer would we obtain if we reversed the orientation of C ? 

[Ans: 1. 27/12; 2. 3/2, while reversing the orientation would give -3/2.] 


3.6 More on line integrals 

View this lesson at http://www.youtube.com/watch?v=S2rT2zK2bdo , [line2]. 



Motivation 





A line integral involves the integration of functions over curves. Applications include: calculating work 
done; determining the total mass and center of mass of thin wires; and also finding flux over curves. 
An important part of the integration process is to appropriately describe the curve of integration by 
^ <c parametrization ? . _ 



Example 


< 



For the vector field 


F = 8e x i + cosh zj — y 2 k 


calculate: 

1) V • F (i.e. div F) and V x F (i.e. curl F); 

2) The line integral of F along the straight line from A(0, 1, 0) to F>(ln(2), 1, 2). 
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— The bigger picture 


For a smooth curve C with parametrization r = r (t), a < t < b 

f fds= [ b f(rmwmdt. 

Jc Ja 

For an oriented smooth curve C with parametrization r = r (t), a < t < b 

/ F • dr — / F • T ds 

Jc Jc 

where T is the unit tangent vector in the direction of motion. 



Exercises 


< 



Let C be the curve parametrized by r (t) = t 2 i + y/t] + v^k, from t = 2 to t = 4. If 
F(x, y, z) := x 2 i — zj + 2xyk then compute 


L 


F • dr. 


v 


[Ans: 49013/3.] 

_ ^ _ ) 


3.7 Fundamental theorem of line integrals 

View this lesson at http://www.youtube.com/watch?v=nH5KCrEuvyA , [FTL1]. 



Motivation 


< 



A line integral involves the integration of functions over curves. In certain cases the value of the line 
integral is independent of the curve between two points. In such a case we may apply a “fundamental 
^theorem of line integrals”. 
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Example 


< 



Consider the vector field 


F := (2 xyz + cosx)i + x 2 z] + ( x 2 y + e z )k. 

1) Show that V x F = curl F = 0. 

2) Find a scalar field 0 such that V0 = grad 0 = F. 

3) Hence or otherwise evaluate J c F • dr where the path C is parametrized by 

r(£) := (sin £, cost, t), from t = 0 to t — 2ir. 

\ __ ) 



/ i 

Maastricht University 


Join the best at 
the Maastricht University 
School of Business and 
Economics! 

gjpj* 

• 33 rd place Financial Times worldwide ranking: MSc 
International Business 
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• 2 nd place: MSc Management of Learning 

• 2 nd place: MSc Economics 

• 2 nd place: MSc Econometrics and Operations Research 

• 2 nd place: MSc Global Supply Chain Management and 

Change 

Sources: Keuzegids Master ranking 2013; Elsevier 'Beste Studies' ranking 2012; 

Financial Times Global Masters in Management ranking 2012 


Maastricht 

University is 
the best specialist 

university in the 
Netherlands 

(Elsevier) 

Master's Open Day: 22 February 2014 

www.mastersopenday.nl | 
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The bigger picture 


• If there is a scalar-valued function 0 such that V0 = F then the line integral 

J F • dr 

is independent of the path C from starting point A to ending point B. Furthermore, 

J F • dr = <f>(B) -<j>{A). 

• If curl F = 0 then a 0 will exist such that V0 = F. 



Exercises 


< 



Consider the vector field 


F := yzi + xz] + xyk. 

^Evaluate f c F • dr where C is any path from (0, 0, 0) to (2, 1, 3). _ [Ans: 3.] 


3.8 Flux in the plane + line integrals 

View this lesson at http://www.youtube.com/watch?v=YFIiTuLNu3Sc , [flux]. 



Motivation 





We show how line integrals can be used to calculate the outward flux (flow rate) of a vector field over 
a closed curve in the plane. Such ideas have important applications to fluid flow. 

V_____ J 



Example 


< 



Let C be the ellipse x 2 /4 + y 2 = 1. 


a) Compute an outward-pointing normal vector n to C . 

b) If F(x, y) 4xi + y] then compute the outward flux over C . 

V_1_ J 
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r 


The bigger picture 


A 


For a closed, smooth curve C with parametrization r = r(t), a < t < b a normal vector to 
C may be produced by calculating (or ). 

For a smooth curve C with parametrization r = r(t), a < t < b and with outward 
pointing unit normal vector n we can compute the outward flux over C via 


F • n ds. 



Exercises 





Let C be the circle with centre (0, 0) and radius 2. 

1. Compute an outward-pointing normal vector n to C . 

2. If F (x,y) := (x + l)i + y] the compute the outward flux over C . 

3. Compute 
V#F dA 



where H is the disc (bounded by C ). Compare your answer with (b). 

[Ans: 1 . If r (t) = (2 cos 2 sin t) , 0 < t < In then n = (cos sin t ); 2. 3. 8n.] 
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4 Double integrals and 
applications 


4.1 How to integrate over rectangles 

View this lesson at http://www.youtube.com/watch?v=My6sdEekbHM , [Dbllntl]. 


r 


Motivation 




Double integrals are a generalisation of the basic single integral seen in high-school. Double integrals 
enable us to work with more complicated problems in higher dimensions and find many engineering 
^applications, for example, in calculating centre of mass and moments of inertia of thin plates. 


Example 


Evaluate 


3 r 3 


/ := / (x 2 — 2 xy + 2 y 3 ) dy dx. 


1 J 2 



Click on the ad to read more 
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The bigger picture 


• When evaluating a given double integral, perform the “inside” integral first and then 
move to the outside integral. 

• Fubini’s Theorem (simple version): If f = f(x,y) is continuous on the rectangle 

R {(x, y) : a < x <6, c < y < d] 


then 


n b nb pd 

f{x, y) dx dy = / / f(x,y) dy dx. 

J a J c 


It is generally a good idea to sketch the region of integration so as to better-understand 
the geometry of the problem. 



4.2 Double integrals over general regions 

View this lesson at http://www.youtube.com/watch?v=9cAVY9niDnI , [DbllntG]. 



Motivation 


< 



We now learn how to integrate over more general two-dimensional regions than just rectangles. 

V_ ) 


Example 


Evaluate 


where 


I := / j ( x 2 y + y s ) dy dx 


Vt := {(x, y) : x 2 + y 2 < 1, x > 0, y > 0}. 


"V 


J 
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r 


The bigger picture 


A 


• To form a integral over more general two-dimensional regions we appropriately describe 
its boundary (or edges). 

• When evaluating a given double integral, perform the “inside” integral first and then 
move to the outside integral. 

• It is generally a good idea to sketch the region of integration so as to better-understand 
the geometry of the problem. 



4.3 How to reverse the order of integration 

View this lesson at http://www.youtube.com/watch?v=zxyx73HAtfQ , [Dbllnt2]. 



Motivation 





The order of integration in double integrals can sometimes be reversed and can lead to a greatly 
simplied (but equivalent) double integral that is easier to evaluate than the original one. This “order 
reversion” technique can be used when performing calculations involving the applications associated 
^with double integrals. 


r 


Example 


Evalute the following integral by reversing the order of integration 



"A 
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r 


The bigger picture 


A 


When reversing the order of integration, the region of integration needs to be suitably redescribed. 
Appropriate diagrams will help with understanding the geometry of the problem and hence will assist 
^with the redescription of the region of integration. 



Exercises 


s 



Evalute the following integral by reversing the order of integration 


/ := 



12 


sin(i/ 2 ) dy dx. 


0 J 3x 


[Ans: (1 - cos 144)/6.] 
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4.4 How to determine area of 2D shapes 

View this lesson at http://www.youtube.com/watch?v=wiFUGkHTmOO , [Area]. 



Motivation 


s 

\ 


Double integrals are a generalisation of the basic single integral seen in high-school. They can be used 
to determine the area of two-dimensional shapes, which can be an important part of, for example, 
determining the centroid. 


Example 


Using double integrals, calculate the area of the region bounded by the curves: 

y = V%; y = x\ y = x/2. 


The bigger picture 


Assume that/is continuous on a two-dimensional region H. We can calculate the area of 
D by evaluating 


dA. 


• If D is not a simple region then it may be possible to deconstruct D into simple parts in 
order to calculate the integral(s) involved. 

• It is generally a good idea to sketch the region of integration so as to better-understand 
the geometry of the problem. 



Exercises 





a) Use double integrals to determine the area of the region bounded by the lines: 

y = x 3 ; y = x 2 . 


b) If you are familiar with double integrals in polar coordinates, then prove that the area of 
a circle with radius R > 0 is nR 2 using double integration. 

[Ans: (a) 1/12 u 2 where u 2 means 'square units”.] 
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4.5 Double integrals in polar co-ordinates 

View this lesson at http://www.youtube.com/watch?v=zknGq9CvZTQ , [Polarl]. 



Motivation 


< 



We know from basic calculus that sometimes a change of variables can greatly simplify very complicated 
integrals. We now explore this idea in the more general setting of double integrals under a specific 
^change of variables known as polar co-ordinates. 


Example 


Consider the double integral 


'■1/V2 ry/l^y 


3 y dx dy. 


1) Sketch of the region of integration. 

2) Evaluate the double integral by applying polar coordinates. 




J 


The bigger picture 


When transfering double integrals from Cartesian co-ordinates to polar co-ordinates, the three 
essential substitutions are: 

x = r cos 6 ; y — r sin 9 ; dA = r dr dO. 

It is generally a good idea to sketch the region of integration so as to better-understand the geometry 
of the problem. 



Exercises 

v_J 


s' 



Describe, in polar co-ordinates, the set of points of the unit disc Cl that has centre (0, 0) and that lie 
in the first quadrant (x > 0, y > 0). Hence, evaluate 


V 


/ := 



x 2 + y 2 dA. 


[Ans: 7r/6]. 


4.6 More on integration & polar co-ordinates 

View this lesson at http://www.youtube.com/watch?v=TOa3hgRT7iA , [Polar2]. 
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Motivation 


< 



We know from basic calculus that sometimes a change of variables can greatly simplify very complicated 
integrals. We now explore this idea in the more general setting of double integrals under a specific 
^change of variables known as polar co-ordinates. 


r 


Example 


Evaluate by using polar co-ordinates 


r2 ^2y-y 2 _ 

/ := / / _ x 1 + y 2 dx dy. 

Jo J-yj2y-y 2 

\ _ 




J 
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The bigger picture 


• When transfering double integrals from Cartesian co-ordinates to polar co-ordinates, the 
three essential substitutions are: 

x — r cos 9 ; y — r sin 9 ; dA — r dr d9. 

• It is generally a good idea to sketch the region of integration so as to better-understand 
the geometry of the problem. 



4.7 Calculation of the centroid 

View this lesson at http://www.youtube.com/watch?v=rIVB-OsNIvY , [Centroid]. 



Motivation 


< 



In engineering, the word centroid means the geometric center of an objects shape. If a physical object 
(for example, a thin plate) has uniform density, then its center of mass is the same as the centroid of 
^its shape. Double integrals can be used to calculate the centroid of simple two-dimensional regions. ^ 



Example 


s' 



Calulate the centroid of the region that lies in the first quadrant of the XY-plane, bounded by the circle 


2,2 2 . r\ 

x + y = a , a>(J. 
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The bigger picture 


The co-ordinates (x, y ) of the centroid of a two-dimensional region Cl can be calculated 
by the use of double integrals via: 


x — 


I fa x dA 

fin ^ ’ 


y = 


IL y d _A 

fin dA 


It is generally a good idea to sketch the region of integration so as to better-understand 
the geometry of the problem. Look for symmetry 



Exercises 



\ 

Calculate the centroid of the region lying in the first quadrant which is bounded by the curves: 


y = 0; y 2 

= 2x ; x + y = 4. 


_ [Ans: (x,y) = (64/35,5/7).] J 

4.8 How to calculate the mass of thin plates 


View this lesson at http://www.youtube.com/watch?v=oty3f7wHY M , [Mass]. 



Motivation 

v_J 


< 



A common and important challenge in engineering is that of calculating the mass of a given object. 
Double integrals can be used to calculate the mass of simple two-dimensional objects, such as thin 
^plates or sheets, if the object has a known density function. 



Example 

V_ ) 


< 



A thin plate occupies the region bounded by the curves: 


y — \J a 2 — x 2 ; x = ±a; y = —a. 

If we associate the density function 

p(x,y) := y + a 

^with each point in the plate, then calculate the total mass of the plate. 
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5 Ordinary differential equations 

5.1 Separable differential equations 

View this lesson at http://www.youtube.com/watch?v=uQuBVa5PPy4 , [SepODE]. 



Motivation 





Ordinary differential equations (ODEs) are used to mathematically model a range of phenomena from 
science and engineering. ODEs naturally arise when making simple physical assumptions about the 
model under consideration. By computing a solution to an ODE we can make precise predictions 
about future and past states of the model of interest. 

Separable ODEs are a special type of ODE and have the form 

dy 


dx 


= f(x)g(y). 



Example 





(a) Solve the separable ODE 

dy x 
dx 2 y 


^ (b) If we have the extra information y(0) = 1 then determine the constant in the solution in (a). 


r 


The bigger picture 


>v 


• To solve a separable ODE simply: separate the variables (via divsion and / or 
multiplication) and then integrate both sides. 

• Sometimes some algebraic manipulation is required to identify a separable ODE. 

• If an initial condition is known, then the arbitrary constant from the general solution can 
be determined. 

v_1_ J 
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Exercises 





Solve the following separable ODEs: 


(a) 


dy_ 

dx 


2x 

3(y + l) 2 



2 xy, y( 0) = 1. 


[Ans:(a) y(x) = —1 + ( x 2 + c) 1 / 3 ;^) y(x) = Ae x2 is the general solution (where A is a constant); 
^and y(x) e x2 is the particular solution yielded from the initial condition.] 


5.2 Linear, first-order differential equations 

View this lesson at http://www.youtube.com/watch?v=hIsLGmrpoPo , [LinODE]. 



Motivation 





Ordinary differential equations (ODEs) are used to mathematically model a range of phenomena from 
science and engineering. ODEs naturally arise when making simple physical assumptions about the 
model under consideration. By computing a solution to an ODE we can make precise predictions 
about future and past states of the model of interest. 

Linear ODEs are a special type of ODE and have the form 

^ + p(x)y = q(x). 

_ J 


Example 


a) Solve the linear ODE 

dy , i 2 
dx x 

for the solution y(x) for v > 0. 


b) If we have the extra information y(l) = 0 then determine the constant in the solution in (a). 
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The bigger picture 


Construct the integrating factor u{x) dx 

Multiply both sides of the ODE by u and collapse the left-hand-side to the derivative of 
y ■ u(x), ie £(y ■ u(x)). 

Integrate both sides and solve for y. 

If an initial condition is known, then the constant from the general solution can be 
determined. 



Exercises 


< 



Solve the following separable ODEs: 


(a) ^ - 2xy = e x \ y{ 0) = 1. 


(b) f x + 2xy 


2x. 


[Ans: (a) y(x) - e xl (x + c) is the general solution and the initial condition yields y(x) - e %1 (x + 1); 
^(b) y{x) = 1 + ce~ x2 .] 
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5.3 Homogeneous, first-order ODEs 

View this lesson at http://www.youtube.com/watch?v=UpLQUGBznE4 , [FOHomogODE]. 



Motivation 


< 



Ordinary differential equations (ODEs) are used to mathematically model a range of phenomena 
from science and engineering. By computing a solution to an ODE we can make precise predictions 
about future and past states of the model of interest. 


A homogenous, first order ODE is of the form 




Example 





Solve the homogeneous, first-order ODE 


/ y-x 

y =- 

y + x 

^for y = y(x) with x > 0. 


The bigger picture 


To solve a first-order, homogeneous ODE: 
• Make the substitution 

y 


u — 


X 


compute the derivative y in terms of u 

obtain a separable ODE in u 

solve for u and then backsubstitute for to obtain y. 
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5.4 2nd-order linear ordinary differential equations 

View this lesson at http://www.youtube.com/watch?v=2symahFSBd8 , [2ndHomog]. 



Motivation 





Ordinary differential equations (ODEs) are used to mathematically model a range of phenomena from 
science and engineering. ODEs naturally arise when making simple physical assumptions about the 
model under consideration. By computing a solution to an ODE we can make precise predictions 
about future and past states of the model of interest. 

We shall discuss “the method of undetermined coefficients” for the case 

ay" + by' + cy = 0 

where a, b and c are given constants. These kinds of problems arise in the study of vibrating systems. 



Example 





Solve the ODEs: 


^(a) y" + 4 y' + 3y = 0; (b) y" + 4 y' + Ay = 0; (c) y" + 2 y' + 5y = 0. 



The bigger picture 





To solve the ODE ay" + by' + cy — 0 

• Solve the corresponding quadratic (characteristic) equation a\ 2 + bX + c = 0 for X. 

• If the solutions X x , \ 2 are real and unequal then y(%) — Ae XlX + Be x<2X . 

• If the solutions are equal, ie \ = \ 2 = X then y(x) — Ae Xx + Bxe Xx . 

• If the solutions X x , X 2 are complex conjugate pairs, ie \ l = a + fii, X 2 = a - jii then 

y(x) = e ax (A cos (3x + B sin f3x). 
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Ordinary differential equations 



Exercises 


< 



Solve the following ODEs: 

(a) y" + 3 y' + 2y = 0; (b) y" + 2y' + y = 0; (c) y" + 4 y' + 5y = 0. 


[Ans: (a) y(x) = Ae~ x + Be ~ 2x ; (b) y(x) = Ae~ x + Bxe~ x ; 

(c) y(x) = e~ 2x (A cos x + B sin x).] 

5.5 Nonhomogeneous differential equations 

View this lesson at http://www.youtube.com/watch?atch?v=haAV- 6R5VU , [2ndInHomog]. 



Motivation 





Ordinary differential equations (ODEs) are used to mathematically model a range of phenomena from 
science and engineering. ODEs naturally arise when making simple physical assumptions about the 
model under consideration. By computing a solution to an ODE we can make precise predictions 
about future and past states of the model of interest. 

We shall discuss the case 

ay" + by' + cy = f(x) 

where a, b and c are given constants and / is a reasonably simple type of function. These kinds of 
problems arise in the study of vibrating systems with F representing an external forcing function. 



Example 





Solve the ODEs: 


^(a) y" — 5 y r + 6y — 2x + 3; (b) y" — 5 y r + 6 y = 12e 5x ; (c) y" — 5 y r + 6 y — 10e 2a; . 


51 


Download free eBooks at bookboon.com 
























Engineering Mathematics: YouTube Workbook 


Ordinary differential equations 



The bigger picture 





To solve the ODE ay" + by' + cy — f ( x ) 

• Solve the corresponding homogeneous problem ay" + by' + cy = 0 for y H . 

• Construct a particular solution y p to the nonhomogeneous ODE, based on the form of 
fix) 

• The general solution is then y = y p + y H . 



Exercises 


s' 



Solve the following ODEs: 

(a) y" + 3 y' + 2y = e x ; (b) y" + 2y' + y = 2x (c) y" + 4 y' + 5y = 3. 


[Ans: (a) y(x) = Ae~ x + Be 2x + e~ x /6; (b) y(x) = Ae~ x + Bxe~ x + 2e~ x - 4; 

(c) y(x) = e~ 2x (A cos x + B sin x) + 3/5.] ^ 
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Ordinary differential equations 


5.6 Variation of constants / parameters 

View this lesson at http://www.youtube.com/watch?v=rhKoNLjOhIg > [VarParam]. 



Motivation 



A 


The method of variation of constants is a very powerful technique for solving 2nd-order, linear ODEs. 
It is useful when the method of undetermined coefficents does not apply. 

We shall discuss the case 

ay" + by' + cy = f(x) 


where a, b and c are given constants and/is a given function. These kinds of problems arise in the 
^ study of vibrating systems. 



Example 





Solve the ODE 


y" + 4 y = cosec 2x. 


The bigger picture 


To solve the ODE ay" + by' + cy — f (x) 

• Formulate the general solution y H to the homogeneous ODE ay" + by' -\- cy — 0, that 
is, yH = Ayi(x) + Byz^x ) with^ and y 2 being linearly independent. 

• Construct a particular solution y to the inhomogeneous ODE ay" + by' + cy = f(x) 
by using the y x and y 2 in 

y r = -*(*) [ A* + mix) f Ax 

J aW (x) J aW (x) 

where W(x) := y^y'^x) - y 2 (x)y[(x). 

• The final (general) solution with have the form y - y H + y p . 



Exercises 

v_ J 





Solve: 

y" + y = secx. 

[Ans: y{x) = A cos x + B sin x + (cos x) In | cos x\ + x sin x] ^ 
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Laplace transforms and applications 


6 Laplace transforms and 
applications 

6.1 Introduction to the Laplace transform 

View this lesson at http://www.youtube.com/watch?v=a7BzFQlItKw , [LT1]. 


r 


Motivation 




Laplace transforms are used in solving initial value problems that involve linear, ordinary differential 
equations with constant coefficients. These types of problems usually arise in modelling of phenomena. 
Laplace transforms offer an advantage over other solution methods to initial value problems as they 
^streamline the process and can easily deal with discontinuous forcing functions. 


Example 


Define and denote the Laplace transform of a function f = f(t) by 

P oo 

£{/(* *)} : = / e ~ st f(t) dt = F(s). 

JO 

1) From the above definition, calculate 


C{e zt }, for all s > —2. 

2) If a is a constant then use (i) to make a conjecture concerning 

C{e~ at }, for all s > —a. 


(7.1.1) 


r 


The bigger picture 




V 


• While it is possible to directly use the definition of the Laplace transform to perform 
calculations, it is standard practice to use a table of Laplace transforms. This can save a 
lot of time! 

• The inverse of the Laplace transform is denoted by . Roughly speaking, it “undoes” 
the Laplace transform and vice-versa. From (ii) we deduce 
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Exercises 


< 



a) From the above definition (7.1.1), calculate 
£{5 1 + 3}, for all s > 0. 


b) If c and d are constants, then use (a) to make a conjecture about what 
C{ct + d}, for all s > 0. 

should be. 

c) Use (b) to calculate 



[Ans: (a) 5/s 2 + 3/s; (b) c/s 2 + d/s; (c) t.] 

V___ J 

6.2 Laplace transforms + the first shifting theorem 

View this lesson at http://www.youtube.com/watch?v=WpUbbsULIds , [LT2]. 
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Laplace transforms and applications 



Motivation 


< 



The first shifting theorem provides a convenient way of calculating the Laplace transform of functions 
that are of the form 


/(* *) := e~ at g(t ) 


where a is a constant and g is a given function. 


Example 


Consider the first shifting theorem 

C{e~ at g{t)} = G(s + a ) where G(s) = £{g(t)}. 

Apply the result to calculate: 

a) £{e 2 V}; 

b) c -i 


(s - 1)* 


r 


The bigger picture 




• For the Laplace transform via the first shifting theorem, calculate G(s) from g(t) and then 
shift it by a units in an appropriate manner. 

• For the inverse Laplace transform via the first shifting theorem, identify a and G(s) from 
G(s + a). Then calculate g(t) from G(s) and multiply by e~ at . 
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6.3 Laplace transforms + the 2nd shifting theorem 

View this lesson at http://www.voutube.com/watch?v=UyoxbrQ3eFU , [LT3]. 



Motivation 





Laplace transforms offer an advantage over other solution methods to initial value problems as they 
streamline the process and can easily deal with discontinuous forcing functions. 

The second shifting theorem gives us a way of computing the Laplace transform of certain types of 
discontinuous functions. 


Example 


Consider the second shifting theorem 

C{u{t — c)g(t — c)} = e~ cs G(s ) where G(s) = C{g{t)} 
and u(t — c) is the Heaviside step function. 

Apply the result to calculate: 

a) C{u{t - l)(t - l) 2 }; 

’ 2e~ 4s ' 


b) £ 


-1 


r 


The bigger picture 




• For the Laplace transform via the second shifting theorem, identify c and g(t) and then 
calculate G(s). 

• For the inverse Laplace transform via the second shifting theorem, identify c and G(s). 
Then calculate g(t) from G(s) and shift it to form g(t - c). 



Exercises 


< 



Apply the second shifting theorem to calculate: 


a) C{u(t - 2)(t - 2) 3 }; 



[Ans: (a) 6e 2s / s 4 ; (b) u(t - 2)(t - 2) 4 /24], ^ 
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6.4 Laplace transforms + differential equations 

View this lesson at http://www.youtube.com/watch?v=SNVfyQ4 avY , [LT4]. 



Motivation 


< 



Laplace transforms offer an advantage over other solution methods to initial value problems as they 
^streamline the process and can easily deal with discontinuous forcing functions. 



Example 


< 



Consider the transform of derivatives identities 


= sY(s) - 2/(0), where F(s) = C{y(t)} 

£{y"(t)} = s 2 Y ( s ) - sy(0) - y\ 0 ), 

Apply the result to solve the initial value problem: 

y" + y = u(t- 1 ), 2 /( 0 ) = 1 , y'(0) = 0 . 

___ J 
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The bigger picture 


When applying Laplace transforms to solve initial value problems: 

1. take the Laplace transform of both sides of the differential equation 

2. apply the transform of derivatives identities 

3. incorporate the initial conditions 

4. solve for Y(s) 

5. use inverse Laplace transforms to produce the solution y(t). 



Exercises 





Apply the transform of derivatives identities to solve: 

y" + y = 2u(t-l), 2 ,( 0 ) = 0, y'{ 0) = 1. 

[Ans: y(t) = sin t + 2 u{t - 1)[1 - cos (t -1)].] 
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LAPLACE TRANSFORMS 

£{/(*)} = Jo°° e ~ st f(t)dt = F(s) 


m 

F(s) 

i 

l/s 

t 

1/s 2 

t m , 

m\ /s m+1 

{v > -1) 

r^+iVs^ 1 

e~ at 

l/(s + 0 .) 

sin bt 

6/(s 2 + b 2 ) 

cos bt 

s/(s 2 + b 2 ) 

sinh bt 

b/(s 2 - b 2 ) 

cosh bt 

s/{s 2 - b 2 ) 

sin bt — bt cos bt 

2 6 3 /(s 2 + 6 2 ) 2 

sin bt + bt cos bt 

26s 2 /(s 2 + b 2 ) 2 

t sin bt 

26s/(s 2 + 6 2 ) 2 

te~ at 

l/(s + a) 2 

u(t — c ) 

e _c */s 

e~ at g(t) 

G(s + a) 

tg{t ) 

-G'(s) 

u(t - c)g(t - c ) 

e _C5 G(s) 

m 

sG(s) - <?(0) 

g"(t) 

s 2 G(s) - sp(0) - tf'(O) 

g"\t) 

s 3 G(s) — s 2 <7(0) — sg^O) — g"{fy 

Jo g(r)dr 

G(s)/s. 


In this table, a and b are any real numbers, c is any non-negative real number, m is any non-negative 
integer, v is any real number larger than -1, Y(y + 1) = J 0 °° e~ x x v dx and u(t) is defined by u(t) = 0, 
for t < 0, u( 0) = 1/2; u(t) = 1 for t > 0. 
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7 Fourier series 


The Fourier series of a function/ = f(t) with period 2 L is defined by 

oo 

Sf(t) := a 0 + 


71=1 L 


nnt Tint 

a n cos —— + b n sin —— 

±j ±J 


where 


a ° := f dt 

J —L 

1 f L Tint 

a n := - y j[t) cos — dt 

i 1 f L r , N . n 7 rt 

b n - J f{t) sin — dt 


and /and/' are piecewise continuous on (-L, L). 

7.1 Introduction to Fourier series 

View this lesson at http://www.youtube.com/watch?v=3Ai3dKVPLcQ , [FS1]. 



Motivation 

v_> 





Fourier series are used in solving differential equations that arise in the study of heat flow and 
vibrations. Fourier series provide a means of approximating discontinuous periodic functions over 
^intervals (rather than just near certain points). 


Example 


Calculate the Fourier series of 


f(t) : = { °> 

with/(f) =/(f + 2n) for all t. 


for 0 < t < 7r; 
for t = 0, ±7t; 
for — 7T < t < 0; 


( 8 . 1 . 1 ) 


61 


Download free eBooks at bookboon.com 























Engineering Mathematics: YouTube Workbook 
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Fourier series 


The bigger picture 


• Integration is required to compute the Fourier coefficients a Q , a n and b n and the 
calculations will simplify when/is odd or even. 

• If/is continuous at a point c then Sf(c ) = /(c). 

• If/has a jump discontinuity at c then Sf(c) = [f(c + ) +/(r)]/2. 



Exercises 





Calculate the Fourier series of 


m ■= 


1, for 0 < £ < 1; 
0, for — 1 < t < 0; 


( 8 . 1 . 2 ) 


with f{t) = f{t + 2) for all t. [Ans: Sf(t) = \ + s i n (2^ — l)7rt. Note that/is neither 

even nor odd, but can you see how to write it in terms of an odd function plus a constant? This can 
help to simplify the calculations!] 


7.2 Odd + even functions + Fourier series 

View this lesson at http://www.youtube.com/watch?v=lBjOw6G3dA , [FSOE]. 



Motivation 





Simplifications occur in the construction of Fourier series when dealing with odd or even functions. 


Example 


Let 


f(x) x, 0 < x < 1. 

Extend/as an odd (even) function f o (/) with period 2 and sketch their graphs on [-3, 3]. Find Sf e 
and Sf Q . Show 

.2 


E l 7T Z 

„„ ( 2 fc - 1) 2 = ¥' 
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r 


The bigger picture 


• Integration is required to compute the Fourier coefficients a Q , a n and b n and the 
calculations will simplify when/is odd or even. 

• Simplifications occur when dealing with odd or even functions. 

• If/is continuous at a point c then Sf(c) =/(c). 

• If/has a jump discontinuity at c then Sf(c) = [f(c + ) +/(c )]/2. 




J 



Exercises 



\ 


Calculate the Fourier series of 

f /_\ , = / 3, for 0 < x < 1, 

J ^ \ —3, for — 1 < x < 0, 

where/(v) =f(x + 2). What is the value of S/(0) and S/(l/2)? 

>oo 12 


[Ans: Sf(x) = EZi sin ( 2fc “ 5 S /(°) = °> S A 1 ' 2 ) = 3 -] 


7.3 More on Fourier series 

View this lesson at http://www.youtube.com/watch?v=toBIWQGPRB4 , [FSR]. 



Motivation 



\ 


Fourier series are used in solving differential equations that arise in the study of heat flow and 
vibrations. Fourier series provide a means of approximating discontinuous periodic functions over 
^intervals (rather than just near certain points). 
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Fourier series 


Example 


Let /(x) = |x| for x E [—7r, tt] with/(x + 27t) = /(x). 

1) Make a sketch of/on the interval [—37T, 3tt] . 

2) Calculate the Fourier series of f 

3) Let 


g(x) = 


1, 0 < X < 7T 

— 1, —7T < X < 0, 


with g(x + 27 t) = g(x). By differentiating the Fourier series in b) ii) (or otherwise) obtain 
the Fourier series of g. 

4) What does the Fourier series of g converge to at x = 7r? 

5) Use the Fourier series for g to obtain a series representation for tt/4. 


SIMPLY CLEVER 


SKODA 



We will turn your CV into 
an opportunity of a lifetime 


Do you like cars? Would you like to be a part of a successful brand? 
We will appreciate and reward both your enthusiasm and talent. 
Send us your CV. You will be surprised where it can take you. 


65 


Send us your CV on 

www.employerforlife.com 



Download free eBooks at bookboon.com 


















Engineering Mathematics: YouTube Workbook 


Fourier series 


r 


The bigger picture 


• Integration is required to compute the Fourier coefficients a Q , a n and b n and the 
calculations will simplify when/is odd or even. 

• If/is continuous at a point c then Sf(c) =/(c). 

• If/has a jump discontinuity at c then Sf(c) = [f(c + ) +/(r)]/2. 

• Simplifications occur when dealing with odd or even functions. 




J 



Exercises 


< 



Calculate the Fourier series of 


f(x) := 3(x + 1), — 2 < x < 2 

where/(v) =f(x + 4). 

[Ans: Sf(x) = 3-12 sin ■] ^ 


7.4 Applications of Fourier series to ODEs 

View this lesson at http://www.youtube.com/watch?v=70KZhL7eVN8 , [FS2]. 



Motivation 

V_ J 


s’ 



Fourier series can be useful in solving differential equations that arise in the study of vibrations. They 
provide a means of approximating discontinuous periodic “forcing” functions over intervals (rather 
^than just near certain points). 
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Fourier series 


Example 


Calculate the solution of the initial value problem 

x" + 50x = f(t), x(0) = 0, x'(0) = 0 


where 


1, for 0 < t < 7r; 

f(t) := ^ 0, for t = 0,7r, 27 t; 

— 1, for 7 t < t < 27r; 


with f(t) =f(t + 2n) for all t. Hint: the Fourier sine series of/is 


(8.4.3) 


E 


—(2 k — 1)7r 


sin(2A; — 1)£. 



Maple has been used to plot the solution. 


The bigger picture 


When applying Fourier series to solve initial value problems: 


• If it is not already given, then calculate the Fourier series of the forcing function. 

• Apply the basic theory of differential equations to calculate y H and y . 
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8 PDEs & separation of variables 

8.1 Deriving the heat equation 

View this lesson at http://www.youtube.com/watch?v=K-CI61wV6IQ , [HeartDerive]. 



Motivation 





In the early 1800s, J. Fourier began a mathematical study of heat. A deeper understanding of heat 
flow had significant applications in science and within industry. A basic version of Fourier’s efforts 
is the problem 


a 2 u xx — Ut\ u{0 , t) = 0 = u(L, i); u{x ) 0) = f(x) 


V 


where: u(x, t ) is the temperature at position x at time t ; a 2 is a constant; and/is a given function. 


J 




MAERSK 


I joined MITAS because 
I wanted real responsibility 


The Graduate Programme 
for Engineers and Geoscientists 

www.discovermitas.com 


Real work 
International opportunities 
Three work placements 




a 


I was a construction 
supervisor in 
the North Sea 
advising and 
helping foremen 
solve problems 



Download free eBooks at bookboon.com 
















Engineering Mathematics: YouTube Workbook 


PDEs & separation of variables 


Example 


Derive the heat equation 


u t — a u x 


( 9 . 1 . 1 ) 


where a 2 k/ pc called the thermal diffusivity (in (length) 2 /time) and k , p and c are positive 
constants that depend on the material of the bar. 


The bigger picture 


Our assumptions and method may be summarized as: 

• The total amount of heat H = H(t) in D (say in calories) is 

rx i 

H(t) — / cpu(x,t) dx. 

J Xn 


Fourier’s law of heat flow says that heat flows from hotter regions to colder regions and 
the flow rate is proportional to u x 

dH 


dt 


= ~KU x (x 0 ,t) ~ (~KU x (Xi,t)). 


Differentiation then gives us two expressions for dH/dt which can be equated to form the 
heat equation. 


8.2 Heat equation & separation of variables 

View this lesson at http://www.youtube.com/watch?v=Ifn6AljZF5g , [SepVar]. 



Motivation 





The method of “separation of variables” is a powerful technique that enables us to solve linear partial 
differential equations that arise in the modelling of various physical processes. In some problems, 
^Fourier series are an important element of the approach. _^ 
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Example 

v_J 


< 



Solve 


4m xx M/ , 

0 < x < 7r, t > 0; 

(9.2.2) 

u( 0 , t) = 0 , 

t > 0 ; 

(9.2.3) 

u( 7T, t) = 0 , 

t > 0 ; 

(9.2.4) 

u(x, 0 ) = 

f(x), 0 < X < 7T. 

(9.2.5) 


with f(x) = sin v - 2 sin 2x. 

\ _ ) 


r 


The bigger picture 


Our method may be summarized as: 




• Assume a solution is in the form u(x, t) = X(x)T(t). 

• Formulate and solve two ODEs for X and T determining which values of the separation 
constant y are relevant. 

• If required incorporate Fourier series. 

V___ J 



Exercises 


< 



Solve (9.2.2)-(9.2.5) but with f(x) = — sin 2x + 3 sin 5x as the initial temperature. [Ans: 

(u(x, t ) = —e~ 16t sin 2x + 3e _200t sin 5x .] 

V___ J 


8.3 Heat equation & Fourier series 

View this lesson at http://www.youtube.com/watch?v=4jYgpHLP6xs , [HeartFourier]. 



Motivation 


< 



The method of “separation of variables” is a powerful technique that enables us to solve linear partial 
differential equations that arise in the modelling of various physical processes. In many problems 
^Fourier series are an important element of the solution method. 
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PDEs & separation of variables 



Example 

v_J 


< 



Solve 


xx U j t ■) 

0 < X < 7T, 

t > 0; 

(9.6) 

u( 0, t) = 0, 

t > 0; 


(9.7) 

u( 7T, t) = 0, 

t> 0; 


(9.8) 

u(x, 0) = x 

VI 

VI 

o 

7r. 

(9.9) 


In previous lectures, we have applied the method of separation of variables to the above problem, 
obtaining 

oo 

u(x, t ) = b n e~ 4n2t sin nx. 

n=l 

where b n are constants to be determined. In particular, (?) yields 

oo 

n(x, 0) = x = b n sin nx, 0 < x < n. 

\ _ 1 1=1 _ ) 


r 


The bigger picture 


Our method may be summarized as: 


“A 


• Assume a solution is in the form u(x, t) = X(x)T(t). 

• Formulate and solve two ODEs for X and T determining which values of the separation 
constant y are relevant from the associated boundary conditions. 

• If required incorporate Fourier series. 

v___ J 


8.4 Wave equation and Fourier series 

View this lesson at http://www.youtube.com/watch?v=4UtGyd86FjI , [Wave]. 
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PDEs & separation of variables 



Motivation 


< 



The method of “separation of variables” is a powerful technique that enables us to solve linear partial 
differential equations that arise in the modelling of various physical processes. In some problems, 
Fourier series are an important element of the approach. We now apply these ideas to solve the wave 
^equation. _ 


Example 


Solve 


9 u xx — u tt , 0 < x < 7r, t > 0; 


u{ 0, t) — 0, u{7r, t) = 0, t > 0; 


u(x, 0) = f(x), U t (x, 0) = 0, 0 < X < 7T. 


(9.4.10) 

(9.4.11) 

(9.4.12) 


where 


f( x ) : = 


x , for 0 < x < 7t/ 2; 
7T — X, for 7 t/2 < X < 7T. 


r 


The bigger picture 


Our method may be summarized as: 




• Assume a solution is in the form u(x, t) = X(x)T(t). 

• Formulate and solve two ODEs for X and T determining which values of y are relevant 
from the associated boundary conditions. 

• If required incorporate Fourier series. 

V___ J 
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